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A diagrammatic technique for two-particle vertex functions is used to describe system- 
atically the influence of spatial quantum coherence and backscattering effects on transport 
properties of noninteracting electrons in a random potential. In analogy with many-body 
theory we construct parquet equations for topologically distinct nonlocal irreducible ver- 
tex functions into which the local one-particle propagator and two-particle vertex of the 
coherent-potential approximation (CPA) enter as input. To complete the two-particle 
parquet equations we use an integral form of the Ward identity and determine the one- 
particle self-energy from the known irreducible vertex. In this way a conserving approxi- 
mation with (Herglotz) analytic averaged Green functions is obtained. We use the limit 
of high spatial dimensions to demonstrate how nonlocal corrections to the d — oo (CPA) 
solution emerge. The general parquet construction is applied to the calculation of vertex 
corrections to the electrical conductivity. With the aid of the high-dimensional asymp- 
totics of the nonlocal irreducible vertex in the electron-hole scattering channel we derive 
a mean-field approximation for the conductivity with vertex corrections. The impact 
of vertex corrections onto the electronic transport is assessed quantitatively within the 
proposed mean-field description on a binary alloy. 



I. INTRODUCTION 

Randomness in the chemical composition or due to defects in solids causes elastic scatterings of charge 
carriers that influence significantly their motion and reduce their mobility. To describe the effects of 
randomly distributed scatterers reliably one has to use approximations taking satisfactorily into account 
quantum coherence between scattering events. In particular, self-consistence between one and two-particle 
functions is needed if we want to assess the role of backscatterings and the backflow on transport properties 
of electrons in a random potential. 

A simplest possibility to account for quantum coherence is to sum all multiple scatterings on a single 
scatterer. A self-consistent theory with all single-sita,scatterings is the Coherent-Potential Approxima- 
tion (CPA) and was developed in the end of sixties.Era It is a mean-field approximation for disordered 
(noninteracting) electrons. It provides a local coherent potential (self-energy) that, in the thermodynamic 
limit, comprises the effects of the random potential on the motion of the single electron. The CPA, as 
other mean-field theories, however, suppresses spatial coherence between distinct scattering centers and 
the moving electrons feci the influence of the random potential only via an averaged medium. 

CPA has proven very successful in the description of one-electron properties at the model level as 
well as in realistic calculations of random alloys.trQ Although this approximation can as well describe 
two-particle averaged functions, due to the lack of spatial coherence, it fails to capture backscattering 
effects on the transport coefficients and the electrical conductivity. The two-particle CPA vertex does 
not depend on the transfer momentum between the incoming and the outgoing particles and hence the 
CPA conductivity reduces in the single-band bulk systems to a contribution from a single particle-hole 
bubble. There are no vertex corrections to the electrical conductivity within the standard CPAB 

Vertex corrections to the single-bubble one-electron conductivity are important in various situations. 
In low dimensions (d < 2) or for sufficiently strong disorder they lead to Anderson localization^ Further 
on, tunnel conductance or transport through multilayered dirty metals are essentially influenced by vertex 
corrections.Q To obtain more realistic results for the electronic transport in dirty metals one has to go 
beyond the standard CPA to the conductivity and develop approximations containing spatial quantum 
coherence and backscattering contributions. 
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There, is a long history of efforts to improve upon the mean- field CPA description of disordered 
electronsu. Most of them concentrate on single-particle properties and improve upon the CPA in the self- 
energy (coherent potential). A natural extension of the single-site theory is to use clusters self-cpjasistently 
embedded in an averaged medium. However, apart from the traveling-cluster approximation|ffiJ exten- 
sions in the lattice^pace fail to warrant global analytic properties of the solution and hence spurious 
effects can emerge.Eail3 Only recently a cluster expansion in momentum space was suggested that war- 
rants [analytic (Herglotz) properties of the resulting averaged propagators and the self-energy at each 
stage. E£l Cluster approximations with self-energy diagrams improve also two-particle vertex functions. 
However, cluster approximations reduce spatial quantum coherence only to a discrete set of lattice sites 
or momenta. Such approximations then remain perturbative in the coherence range and cannot lead 
to Anderson localizMion to which we need long-range coherence with infinite-many backscattering or 
"crossed" diagrams OuJ 

Using cluster approximations to improve upon the mean-field transport properties means that we have 
first to extend the one-electron calculation scheme. A tremendous effort at the one-particle level is to 
be exerted to obtain significant changes in transport properties. Cluster expansions are hence not very 
effective in calculating quantum coherence effects in the electrical conductivity. It is more efficient to 
develop approximations directly for two-particle functions. 

A suitable framework for developing two-particle approximations is a parquet approach devised within 
quantum many-body theory.Ejtj It is an advanced scheme of summation of Feynman (many-body) 
diagrams based on renormalizations of two-particle vertex functions. Its main idea is to utilize ambiguity 
in the definition of the two-particle irreducibility. Each two-particle irreducibility, i. e., the way how 
pairs of one-particle propagators are cut without disconnecting the diagram, defines a scattering channel 
and a Bethe-Salpeter equation for the full two-particle vertex. Since different two-particle channels are 
topologically incquivalent, a solution of the Bethe-Salpeter equation from one channel (reducible function) 
is irreducible in the other channels where it is used in the input (integral kernel) in the respective Bethe- 
Salpeter equations. Thereby a set of coupled, nonlinear self-consistent equations for the two-particle 
irreducible vertices (parquet equations) is obtained. Parquet equations have been applied onto various 
many-body problems, but no significant attempt has been made to use the parquet-type renormalization 
of Feynman diagrams in disordered systems. 

In this paper we develop a parquet approach to systems with noninteracting electrons subjected to 
a random potential. We show how to construct controlled approximations directly for the two-particle 
vertex using the idea of parquet diagrams. Since the parquet construction applies only to nonlocal propa- 
gators, we start from the limit of high spatial dimensions where the diagonal and off-diagonal one-particle 
propagators separate and the CPA becomes exacted Beyond this limit we construct parquet equations 
for two-particle irreducible vertices from Bethe-Salpeter equations with a perturbed nonlocal one-particle 
propagator and the local two-particle vertex as input. Next we use a Ward identity to determine the 
self-energy and the full one-particle averaged propagator from the calculated vertex functions. This 
self-consistent procedure warrants conservation laws and analytic properties of the one-particle functions 
whenever the solutions to the two-particle parquet equations are analytic. 

The unrestricted system of parquet equations is not soluble in general. We hence resort to high 
spatial dimensions where the two-particle self-consistence-js naturally suppressed and one obtains the 
asymptotic form of the two-particle vertex in closed formO We use this explicit result to derive a mean- 
field approximation for the electrical conductivity containing vertex corrections. We then choose a binary 
alloy to make quantitative assessments of the impact of vertex corrections on the bulk conductivity. 
The paper is organized as follows. We derive in Sec. [n] the parquet equations for vertex functions of 



disordered electrons. In Sec. Ill we show how Ward identities can be used to determine the self-energy 
from a given irreducible vertex function so that we preserve conserving character of the approximation. 
The electrical conductivity with the irreducible vertex function in the electron-hole scattering channel is 



derived in Sec. IV where we use the result to constructing a mean-field approximation for the electrical 
conductivity with vertex corrections. The mean-field approximation is applied in high spatial dimensions 
on a binary alloy in Sec. |v| to obtain quantitative results. 

II. CALCULATION OF THE VERTEX FUNCTION 

We use a tight-binding Anderson model of noninteracting spinless electrons moving in a random, 
site-diagonal potential Vi described by a Hamiltonian 
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Had = E ^' c l c j + E V ^ c * 



(1) 



The values of the random potential Vi are site-independent and obey a disorder distribution p(V). I. e., 
a function depending on the random potential Vi is averaged via 



(X(Vi)) av = / dVp(V)X(V) . 



(2) 



Solving the problem of disordered electrons in thermodynamic equilibrium amounts to finding the 
averaged free energy defined as 



F av = -k B T (lnTr exp^pHADitijVi)} 



(3) 



where the trace Tr runs over the electronic degrees of freedom in the Fock space. However, the aver- 
aged free energy does not contain the entire information about the disordered system. In particular we 
cannot derive transport properties and the response to disturbing external forces from it. We need to 
know averaged products of Green functions for different energies. To include external perturbations into 
the thermodynamic description we introduce a new quantity Q v (Ei, E 2 , ■ ■ ■ , E v \ U). It is a generalized 
averaged grand potential with v energy states coupled via an external perturbation U. We define 



n v {E 1 ,E 2 ,...,E„;U) = -k B T(lnTrexp{ -(3 J2 [ H 

'•j 1 



W X.. — R. mO)*. ■ 



AH^ 



(4) 



where we assigned to each (complex) energy Ei a separate Hilbcrt state space and AH^' 



Potential ^(Ei, E2, . ■ ■ , E v ; U) is a generating functional for averaged products of 
Green functions up to the z^th order. In practice, within linear-response theory we will use only one and 
two-particle Green functions, i. e., fl v (E±, E2, ...,E V ;U) is expanded up to U 2 . 

Averaged Green functions (propagators) are fundamental quantities with the aid of which we can 
calculate all characteristics of the disordered system. We can use momenta as good quantum numbers, 
since translational invariance is restored for the averaged quantities. Averaged propagators can then be 
expressedpas sums of Feynman diagrams for disordered systems in analogy to the standard many-body 
diagrams .E3 

The averaged one-particle propagator is represented with the aid of the the self-energy or a coherent 
potential E(k, z) that comprises the influence of fluctuations of the random potential onto the motion of 
the single electron. We write 



G(k,z) 



-e(k)-E(k,z) 



1 

N 



-ik(R 4 — Rj) 



zl-t-V 



(5) 



where the first equality expresses the Dyson equation connecting the irreducible one-particle function 
(self-energy) with the one-particle averaged propagator. 
The averaged two-particle propagator is defined as 



G (2) 



ij.kl 



(zi,z 2 ) = 



Zxl-t-V 



z 2 l-t-V 



to which we define the Fourier transform to momentum space as follows 

(k ls Zl , k a , z 2 ; q) = ^E e-^e^^ e-^+^^^Gf^, z 2 ) . 



(G) 



(7) 



ijkl 



The two-particle Green function G^ 2 -*(ki, z%, k. 2 , z 2 \ q) contains also uncorrelated motion of two separate 
particles. The actual measure of a correlated motion of two particles is a vertex function 



T (ki , zi , k 2 , z 2 ; q) - G- 1 (ki , zi ) G~ 1 (k 2 , z 2 ) [g< 2 » (ki , z x , k 2 , z 2 ; q) 

-<y(q)G(ki, Zl )G(k 2 ,z 2 )} G-^ki + q, Z 1 )G- 1 (k 2 + q, z 2 ) . 



(8) 
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In analogy to the Dyson equation for the one-particle propagator we can represent the two-particle 
vertex with the aid of a two-particle irreducible vertex A and a Bethe-Salpeter equation. Unlike the one- 
particle case, the Bethe-Salpeter equation is not defined unambiguously whenever we work with nonlocal 
propagators. This fact we utilize later in the parquet construction. For the present moment we take 
the Bethe-Salpeter equation in the electron-hole channel describing multiple scatterings of a pair of an 
electron and a hole and write 



r(ki,2i,k 2 ,22;q) = A(ki, zi, k 2 , z 2 ; q) + A(ki, z x , k 2 , z 2 ; q") 

q" 

xG(ki + q", Zl ))G(k 2 + q", z 2 )r(k! + q", Zl ,k 2 + q", z 2 ; q - q") . 



(9) 



The one and two-particle irreducible functions, i. e., self-energy £(k, z) and vertex A(ki, z%, k 2 , z 2 ; q) 
are the quantities that we have to approximate in order to determine the one and two-particle charac- 
teristics of a disordered system. The two functions are not completely independent. In a conserving and 
thermodynamically consistent approximation we have a generalized differential Ward identitycd 



A ki, zx, k 2 , z 2 ; q = — — ■ ■ — 

dG(ki + q, zx, k 2 + q, z 2 ; U) 



(10) 



u=o 



We could use it for the determination of the irreducible vertex if we knew explicitly the self-energy as 
a functional of the averaged propagator in the presence of the external disturbance U . It is rarely the 
case. We, however, show in Sec. [II how to use an integral form of the Ward identity to determine the 
self-energy from the known irreducible vertex A. It is then sufficient to construct an approximation for 
the two-particle irreducible vertex A, which will be done in the following subsections. 



A. Local approximation 



We start building approximations to the two-particle vertex function from a local solution where we 
completely loose momentum dependence. The local approximation means that we use only site-diagonal 
one-particle propagators in the perturbation diagrammatic expansion for the irreducible functions. The 
local approximation is best derived within the, Baym-Kadanoff rcnormalized perturbation expansion in 
the limit of high spatial dimensions d — > oo.E!l In this limit we have the following asymptotics for the 
one-particle functions 



G = G dia9 [d } + G o// [fr 1/2 ] , E = Y, dia9 [d°] + E o// [d 



-3/2i 



(11) 



leading to separation of the diagonal (local) and off-diagonal (nonlocal) parts. In the strict limit d = oo 
we can completely neglect the off-diagonal elements and recover the CPA for the self-energy. The defining 
equation in the presence of the external local disturbance U reads 



G{z\,zi\ U) 



G- 1 (z 1 ,z 2 ;U) + Z(zi,z 2 ;U)-Vi 



(12) 



where G(z 1; z 2 ; U) = N~ 2 X)kik 2 ^0*1) z i> z 2] U) is the local element of the matrix one-particle propa- 
gator. The matrix character is forced by the external disturbance which mixes different complex energies. 
Since we are interested only in averaged two-particle functions in equilibrium, we can resort to two energies 
and a two-by-two matrix 



G- i (k 1 ,z 1 ,k 2 ,z 2 ;C/) = 



_(zi- e(ki) - £(zi) U - E(zi, z 2 ; U) 



U -T.(z2,zi\U) z 2 - e(k 2 ) - E(z 2 ) 



(13) 



where e(k) is the lattice dispersion relation. 

We use the Ward identity ( |Io| ) to determine the two-particle irreducible vertex in equilibrium 



A(zi, z 2 ) 



6Y,u(zi, z 2 ) 



8Gu(zi,z 2 ) 
1 

G( Zl )G(z 2 ) 



[7=0 



1 - 



1 + (E(*x) - Vi) G( Zl ) 1 + (E(z 2 ) - Vt) G{z 2 ) 



(14) 
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The full vertex function is then determined from the Bethe-Salpeter equation (|9|) where the one-particle 
propagators are replaced with the local ones. The integral equation reduces to an algebraic one and we 
obtain an explicit representation 

i \ A(z 1; z 2 ) . . 

1 - A(2i, z 2 )G(zi)G(z2) 

For later convenience we denoted the full local two-particle vertex with the lower-case 7. 

Note that the above local approximation coincides with the CPA only for the self-energy and the 
local part of the vertex function. The CPA two-particle vertex is obtained from Eq. (^|) where only the 
irreducible vertex A(zi, z 2 ) is assumed local, but not the one-particle propagators. □ Here we deliberately 
separated the contribution from the off-diagonal propagator. We describe it within the parquet approach 
as a correction to the local approximation to the vertex function. 



B. Nonlocal contributions: Parquet approach 



To go beyond the local approximation we have to distinguish two one-particle propagators. First, 
we have the propagator from the local approximation that we denote G ioc (k, z), from which we actually 
need only the local element G loc (z) = N^ 1 ^ k G io f(k, z). This propagator is defined by the Dyson 
equation with the local self-energy £(z) from Eq. ( pq ) with [7 = 0. Second, we have to introduce a 
new propagator G(k, z) that is defined by the Dyson equation with a nonlocal self-energy E(k, z) to be 
determined later. It is treated in the equations for the two-particle vertex as an external function with 
appropriate analytic properties. In the expansion beyond the local approximation we use a perturbed 
propagator G(k, z) = G(k, z) — G loc (z) to avoid multiple summations of single-site diagrams contained 
in the local approximation. 

We classify nonlocal contributions to the two-particle vertex by the type of the correlated two-particle 
propagation. We either simultaneously propagate an electron and a hole or two electrons (holes). Dia- 
grammatically it means that we connect spatially distinct two-particle scattering events with antiparallel 
or parallel pairs of one-particle propagators. Multiple scatterings of pairs of the same type define a channel 
of two-particle irreducibility. We call a diagram two-particle irreducible if it cannot be split into separate 
parts by cutting simultaneously either electron-hole or electron-electron propagators. The two definitions 
of the two-particle irreducibility lead to topologically inequivalent irreducible functions and to different 
Bcthc-Salpctcr equations for the full vertex. In each Bcthc-Salpctcr equation the two-particle functions 
are interconnected via one-particle propagators in a different manner. We can generically represent the 
Bethe-Salpeter equations as 



r(ki,2i,k 2 ,z 2 ;q) = A Q (ki, z x , k 2 , z 2 ; q) 



A a GG r 



(ki,zi,k 2 ,z 2 ;q) 



(16) 



We used © for the channel-dependent multiplication of the two-particle functions. Here A" is the irre- 
ducible vertex in the a-channel. Although the irreducible functions are different in different channels the 
solution of the Bethe-Salpeter equations must always be the same, the full two-particle vertex T. The ma- 
trix multiplication in momentum space in the electron-hole and electron-electron channels, respectively 
explicitly reads 



XGG»Y 



(ki,2i,k 2 ,2 2 ;q) = — ^X(ki,2i,k 2 ,£ 2 ;q") 



N 



q" 



xG(ki + q", Zi)G(k 2 + q", z 2 )Y(k 1 + q", z 1; k 2 + q", z 2 ; q - q" 



XGGoY 



(ki,zi,k 2 ,z 2 ;q) = — ^X(k x , zi, k 2 + q", z 2 \ q- q") 



(17a) 



xG(kj + q q", Z!)G(k 2 + q", z 2 )Y{k 1 + q - q", z u k 2 , z 2 ; q") . 



(17b) 



Electron-hole and electron-electron channels are not the only inequivalent representations of multiple 
two-particle scatterings. If we allow for hopping between distant sites we no longer can distinguish 
between multiple scatterings of distinct electron-hole pairs or nonlocal self-scatterings of a single particle, 
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i. e., scatterings between the incoming and the outgoing particle at the two-particle vertex. We hence 
have to introduce a third two-particle irreducibility and Bethe-Salpeter equation. We call this third type 
of two-particle scatterings "vertical channel" , since we renormalize only one line of the pair and the ladder 
of multiple scatterings grows "vertically" above (below) the two-particle vertex. We can write the third 
multiplication scheme for two-particle quantities as follows 



XGG + Y 



1 



(ki,zi,k 2 ,z 2 ;q) = — ^ X(ki, z u ki + q, zi; k" - k x ) 

k" 

xG(k", Zl )G(k" + q, zi)F(k", z u k 2 , z 2 ; q). 



(17c) 



In this case only the particle with energy z\ is renormalizcd. Analogously we introduce self-scattering 
vertex corrections to the particle with energy z 2 . The vertical channel actually splits into two, upper and 
lower parts according to whether we renormalize the first or the second energy in the pair. 

Using the above multiplication schemes for different two-particle channels we can write down the 
corresponding Bethe-Salpeter equations explicitly. We choose a subscript a = ± for complex half-planes 
from which we take the particle energy. We use the standard diagrammatic representation for these 
equations and obtain for the electron-hole and electron-electron channels 



ak 



ak + q 




a'k' 



\ / 

a'k' + q 




+ 




k' + q" 



(18a) 



ak ak + q 

X / X / 



k + q - q" 



Aee 
eta' 



+ 



a'k' 



a'k' + q 




k' + q" 



(18b) 



The Bethe-Salpeter equation in the vertical channel is split into two. First we account for self-scattering 
vertex corrections to the upper line and then the same for the lower one. We then have 

ak 



ak 



a'k' 



ak + q 



a'k' + q 



A" 



+ ak + q 



k" + q 




a'k' 



a'k' + q 



(18c) 



G 



ak 



a'k' 



a'k' + q 



ak + q 



ak + q 

/ \ 




a'k' 



a'k' + q 



(18d) 



Note that it is the irreducible vertex from the electron-hole channel A that determines the kernel of the 



integral equations (18c) and (18d). The new vertex A" enters only via the absolute term. The irreducible 
vertices A" and A eh are not generally equal and it is useful for our construction to distinguish the vertical 
channel from the other two. 

In all these equations the double- prime momenta are summed over. The full lines stand for the 
perturbed propagators G(k, z). Otherwise we would encounter multiple summations of single-site scat- 
terings. Separation of the local and nonlocal contributions is important, since in the local approximation 
the three channels coincide. This can be seen from Eq. ( fl7j ) when we insert a local propagator. This fact 
is physically obvious, since we cannot distinguish between the particle and the hole in single-site multiple 
scatterings. The electron and the hole propagators equal. There is only one two-particle irreducible 
vertex A(zi, z 2 ) in the local approximation (CPA). 

We now use the topological inequivalence of the three channels. This means that a reducible function 
from one channel is irreducible in the other ones. If we denote I the completely irreducible two-particle 
vertexCJ we can write 



A a (ki,zi,k 2 ,z 2 ;q) = i"(ki,zi,k 2 ,2 2 ;q) 



E 



(r< 



ki,zi,k 2 ,z 2 ;q) - A" (ki, z\, k 2 , z 2 ; q) 



(19a) 



since the reducible function is a difference between the full and the irreducible vertex from a given channel. 
We use the Bethe-Salpeter equations (|l8|) on the r.h.s. of Eq. (19a) from the corresponding a' channel to 



get rid of the full vertex. We then obtain a set of equations for the channel-dependent irreducible vertices 
A". The completely irreducible vertex is input to these so-called parquet equations. The lowest-order 
contribution to the input function is the local two-particle vertex 7. The parquet equations in a generic 
form are 



A Q (ki,zi,k 2 ,z 2 ;q) = ~/(z 1 ,z 2 ) + ^ j 

a' 



1 - 



A Q GG 



of 



A Q GG © A c 



(ki,zi,k 2 ,z 2 ;q) . (19b) 



Special attention should be pa id t o the v ertic al channel, a! = v, where we have a two-st ep c onstruction 
and Bethe-Salpeter equation (18c) and (18d). The second term on the r.h.s. of Eq. ( |l9b| ) should be 

replaced by jl - A eh GG *} A"jl-* GGA eh X -A". 

Equations ( |l9| ) constitute the parquet approximation for the irreducible vertices for given one-particle 
propagators G(k, z), G loc (z) and the local two-particle vertex 7(2:1, z 2 ), that again is a function of G loc 
from Eq. @. The latter propagator is determined from the local approximation, CPA. The former 
is treated here as a function of a self-energy E(k, z) which we connect to the vertex functions later 
on. Parquet equations ( |19b[ ) form a set of nonlinear integral equations self-consistently determining the 
two-particle irreducible vertex functions A a [S; G loc ]. 

The parquet equations represent a substantial extension of the local approximation. They are, how- 
ever, much more difficult to solve than cluster or other short-range extensions of the CPA at the one- 
particle level. To get a feeling how a solution to the parquet equations looks li ke w e resort to the 
asymptotic limit d —> 00 where the corrections to the local vertex 7(2:1, z 2 ) in Eq. (19b) asymptotically 
vanish.EU The two-particle self-consistence becomes asymptotically insignificant and such a situation can 
be dealt with exactly. 



7 



C. Asymptotic limit of high spatial dimensions 



We know that the off-diagonal, nonlocaLelements are scaled to zero in high dimensions and loose 
their weight with respect to the local ones.c3 However, when summed over the lattice sites they can 
produce finite contributions even in the limit d — oo. It is the case of the two-particle vertex. In 
the leading asymptotic order the irreducible vertices in the Bethe-Salpeter equations become local and 
coincide with 7(21,22). The asymptotic behavior of the full vertex then depends upon which matrix 
element we calculate. We obtain different asymptotic solutions in different channels that we denote 
pa|2-l The corresponding asymptotic Bethe-Salpeter equations in the three channels have the following 
diagrammatic representation 



ak 



ak + q 



k + q" 





+ 



/ \ / \ X 

a'k' a'k' + q 




ak 



ak + q 



pee 
aof 



a'k' 



+ 



a'k' + q 



ak 



k' + q" 



k + q - q" 



k' + q" 



(20a) 




(20b) 



ak 



a'k' 



ak + q 



a'k' + q 



+ ak + q 



k" + q 




a'k' 



a'k' + q 



(20c) 



ak 



ak + q 



ak 



ak + q 



aa' 



a'k' 



a'k' 




a'k' 



a'k' + q 



(20d) 



The one-particle propagator in the high-dimensional Bethe-Salpeter equations equals the CPA solution, 
G(k,z) = G /oc (k, z), or E(k, z) — TP PA (z). We sum over the intermediate momenta in the Bethe- 
Salpeter equations and also the nonlocal part of G loc (k, z) contributes to the leading asymptotic behavior 
of the solution for fixed two-particle momenta. 

The Bethe-Salpeter equations (|2Cj ) in high dimensions become algebraic and can be solved in closed 
form. If we denote 



y ± (q; Zl ,z 2 ) = G off (k, zi)G o// (q ± k, z 2 ) 

k 

the solutions in the three channels are 
r e ' l (q;zi,z 2 ) 

r ee (q;zi,z 2 ) 



l{zi,Z2) 

1 - y(zi, z 2 )Y+(q; z x ,z 2 ) 



l-~/(zi,Z 2 )Y (q;z 1 ,z 2 

2 



r°(q; zuz 2 ) = l(zi,z 2 ) J| 



1 



1 [l-liz^ZijY+^z^Zi)] 



(21) 

(22a) 
(22b) 

(22c) 



The three different vertex functions represent the three asymptotic solutions we obtain in the limit of 
high spatial dimensions. Each solution is characterized by a two-particle momentum that is kept fixed. 
They are in the notation from the preceding subsection ki — \t 2 , ki + k2 + q, and q for the electron- hole, 
electron-electron, and the vertical channel, respectively. The full vertex T in high dimensions reduces to 
a sum 

r(ki, zi, k 2 , z 2 ; q) = T e/l (ki - k 2 ; zi, z 2 ) + T ee (ki + k 2 + q; zi,z 2 ) + r"(q; z u z 2 ) - 2 1 {z u z 2 ) (23) 

where we had to subtract appropriately the local vertex to avoid double counting. Note that the CPA 
vertex as derived by VelickyO equals T eh and does not correspond to the leading high-dimensional asymp- 
totics of the exact two-particle vertex. 

Representation (|2^) uses the natural high-dimensional separation of the diagonal and off-diagonal 
parts of the one-particle propagator. We can rewrite the asymptotic solution (^2|) to another form with 
full one-particle propagators, which is more appropriate for finite-dimensional systems. We denote a 
two-particle bubble 



± (q; z\, z 2 ) = — ^2 G ( k > z i) G (q ± k - z 2) ■ 



(24) 



If we further use relation we obtain a new representation for the two-particle vertex 
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r(ki,zi,k 2 ,2 2 ;q) = A(zi,z 2 ) 



1 



1 - A(zi, z 2 )x+(k 2 - ki; z u z 2 ) 1 - A(z x , z 2 )G(z 1 )G(z 2 ) 



+ 



1 - A(zj, Zi)G(zi)G(zi) 



n 



i 



1 - A(zi, z 2 )G(z 1 )G(z 2 ) 1 - A(zi,z 2 )x~(ki +k 2 -q;zi,z 2 ) 



(25) 



where A(z\, z 2 ) is the CPA irreducible vertex. 

Note that the asymptotic solution for the two-particle vertex can no longer be represented via a single 
Bethe-Salpeter equation. It is because we have three asymptotically equally important but topologically 
different contributions. Each of them is marked by a different fixed two-particle momentum. 



III. WARD IDENTITIES AND THE SELF-ENERGY 



In the preceding section we derived a closed set of parquet equations determining the irreducible two- 
particle vertex functions from the given propagators G loc (z) and G(k, z). The former is the site-diagonal 
part of the CPA propagator but the latter has not yet been specified. It is determined by a self-energy 
E(k, z) from the Dyson equation (^J). We have not yet demanded any relation between the self-energy and 
the vertex functions. We, however, know that thermodynamic consistence demands that the irreducible 
one and two-particle functions be not independent. They are related by the differential Ward identity 
(|To|). To turn the parquet equations a conserving approximation we have to fulfill the Ward identity 
and relate the self-energy to the solution of the parquet equations. The functional differential identity 
( |l0| ) is of little practical help in calculating the self-energy from the irreducible two-particle functions. 
Fortunately there are integral forms of the Ward identity relating one and two-particle averaged functions 
that can be used to determine the self-energy from a known two-particle irreducible vertex. 

A first integral Ward identity was derived by Velickyt] in the framework of the coherent-potential 
approximation. It holds quite generally beyond the CPA and reads in our notation 

E G %jk (*i . *a) = - ^ AG lfc , (26a) 

3 

where Az = z\ — z 2 , AG^ = Gij(zi) — Gij(z 2 ). We show in Appendix A that this identity is a consequence 
of completeness of the eigenstates of the Hamiltonian and hence of probability conservation. This identity 
can be rewritten in momentum space 

^E G(2) ( k ' z i' k ' z 2;q)--^[G(k,Z!)-G(k,z 2 )] . (26b) 
q 

Ward identity (|2^) is inconvenient for application, since it connects one and two-particle averaged 
functions. We better had a relation between irreducible one and two-particle functions alike Eq. (|l^). An 
integral form of such a Ward identity was proven diagrammatically by Vollhardt and WolfleEj and reads 

q 

The diagrammatic derivation of the Vollhardt- Wolfle identity (|27]) utilizes a symmetry of the Anderson 
disordered model but does not establish a direct relation to conservation laws as the derivation of the 
Velicky identity ( |26|). We show in Appendix A that the Vollhardt- Wolfle identity with ki = k 2 follows 
directly from Eq. (|26|) and the Bethe-Salpeter equation. At least this simplified form of the Vollhardt- 
Wolfle identity can be shown to be a consequence of probability conservation^ 

To complete the parquet approximation for the two-particle irreducible vertex functions we can use 
Eq. ( p7|) for ki = k 2 and limiting values of the complex energies along the real axis. We chose z\ = lu+ 
and z 2 = u- where lu± = uj ± ir\ and r\ \ 0. We obtain a relation 
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Im£(k,cj + ) = — ^A e,l (k,w + ,k,w_;k-k')ImG(k',cj + ) (28a) 

k' 

determining the imaginary part of the self-energy along the real axis from the known irreducible vertex in 
the electron-hole channel. A similar formula can be derived for the irreducible vertex from the electron- 
electron channel, cf. Appendix A. 

Next we rely on analytic properties of the self-energy in the upper and lower complex half-planes and 
determine its real part along the real axis from the Kramers-Kronig relation 

ReE(k,w+) = P/ , + . 28b 

The self-energy from the cut along the real axis can be analytically continued to energies in the u ppe r 
and lower complex half-planes. The one-particle perturbed propagator from the parquet equations ( pL9 q ) 
can now be defined from the irreducible vertex via 

G(k, u+) = [u+ - e(k) - E(k, UJ+)}- 1 - G 1oc (uj+) (29) 

with the self-energy from Eq. ( ]2S| ) and the local propagator from Eq. (|l2|). 

Equations (p8|), (|29| ) complete the parquet approach to disordered systems and make it a consistent 
and conserving scheme approximating s imul taneously both the one and two-particle irreducible functions. 
Neither the Bethe-Salpeter equations ( |l9bj ) nor Eqs. (|2|) violate analytic properties, unless there is a 
transition to another phase such as Anderson localization. The analyticity of the parquet approximation 
in the diffusive regime is then completely determined— by its input, the two-particle local CPA vertex 
7(zi, Z2 ) that is known to possess Herglotz analyticity .E3 

The advantage of the parquet approach to disordered systems is that we do not need to bother about 
the diagrams contributing to the self-energy so that-tp obtain a consistent approximation with the proper 
analytic behavior of the averaged Green functio ns JHj An explicit sufficient condition for analyticity of a 
solution of the parquet approximation from Eq. ( |28a|) is 

A e,I (k,w + ,k,w_;q) > 0. (30) 

It can be checked in each step of iterations toward the full solution of the parquet approximation for the 
one and two-particle irreducible functions. 



IV. ELECTRICAL CONDUCTIVITY 



The CPA constitutes a rather good approximation for the one-particle self-energy but completely fails 
to incorporate coherence in the propagation of pairs of particles. The emphasis in the parquet approach 
is hence laid on a systematic construction of diagrammatic approximations for nonlocal two-particle 
functions. In particular those determining transport properties and reflecting Anderson localization. In 
this section we show how a solution of the parquet approximation can be used in the calculation of the 
electrical conductivity. 

We use a Kubo formula for the electrical conductivity with the current-current correlation function. 
If o a f3 is the complex conductivity and Ii a p the current-current correlation function we can writeEl 

a Ql3 (q,E + ) = ~n Qf) (t l ,E + ) (31) 

where again E+ = E + i0 + . The current-current correlation function can be expressed via a Kubo formula 
with the full two-particle vertex T 

n Q/3 (q, iu m ) = P° ( £ ( k + «») " £ (- k )) 9 ( £ ( k ' + «0 " £ (- k '))] 

111 k,k' 

oo 

xk B T {G(k,iw„)G(k + q,^„ + w m )[<5(k-k') 

71— — OO 

+ T (k, ioj n ,k + q,i(uj n + v m );\s! - k)) G(k', ioj n )G(k' + q.,iw n + iv m )]} ■ (32) 
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Here Lo n = (2n + l)irT and v m = 2rmrT are Matsubara frequencies. 

We are actually interested only in the real part of the complex conductivity for real energies. We can 
then analytically continue the expression on the r.h.s. of Eq. ((32|). For the real part of the conductivity 
we obtain 



Re a af} (q, E) = - -jL ^ (v a (k + q) + v a (k) ) (v a (k' + q) + v (k') ) 1 ^ a ^ ^ r 

k,k' cr t 

r ^ It+^tJM ( G(k) w + J(T0 +)G(k + q, ^ + B + zrO+) [*(k - k') 

+r(k, w + iaO + , k + q, w + £ + ?:t0+; k' - k)G(k', lu + iaO + )G(k', lu + irO + )] } (33) 

where f(E) is the Fermi function and we denoted the group velocity v a (k) = m _1 9e(k)/9fc Q and a, r = 
±1. In most situations the static optical conductivity (q = 0, E = 0) at zero temperature is of interest. 
Expression ( |33"|) reduces in this case to 

Re = |^^2 E « Q (k)^(k') ^(-ar)G(k, B F + iaO+)G(k, £ F + 4 r0+) [5(k - k') 

k,k' or 

+ T(k, E F + ia0 + , k, B F + ir0+; k' - k)G(k', B F + iaO + )G(k', E F + ir0+)] . (34) 

We immediately see that the vertex function T contributes to the conductivity only if it depends on 
the transfer momentum k' — k. It is due to the symmetry e(k) = e(— k) and v a (k) — — v a {— k). Since 
the CPA vertex does not depend on the transfer momentum, the vertex corrections to the single-bubble 
electrical conductivity vanish in the CPA treatment. 

Eq. (34) is not appropriate for approximate calculations of the electrical conductivity. The vertex 
corrections contained in T are added to the single-bubble conductivity so that their negative contributions 
may reverse positive sign of the conductivity, thus leading to unphysical behavkw;. 

To avoid such a situation we represent the conductivity in a different way.Eil We utilize the Bethe- 
Salpeter equation in the electron-hole channel to represent the vertex T via the irreducible one, A eh . 
Inserting its formal solution into Eq. ( |34| ) we obtain a new, equivalent representation for the conductivity 

Rea a = ^-^^v a (k)v (k')Y,(-aT)G a (k)G T (k){l- [Af T G CT G r ] (k, k; k' - k) . (35) 



k,k' 



where we denoted G CT (k) = G(k, E F + ia0+) and Af T = T(k, E F + ia0+ , k', E F + ir0+; q). At least for 
not too strong disorder the norm of the operator A^.G CT G T < 1. Hence the conductivity remains in 



this representation non-negative and free of spurious unphysical behavior. However, unlike formula (34), 
representation (|35j) is implicit and its application is conditioned by our ability to solve the Bcthe-Salpeter 
integral equation in the electron-hole channel explicitly. 



A. Asymptotic expression in high dimensions 

To be explicit in the assessment of the contributions of the two-particle vertex to the electrical conduc- 
tivity we again resort to the limit of high spatial dimensions where we know the vertex and the self-energy 
explicitly. We use a hypercubic lattice with a dispersion relation e(k) = — i/(2d) 1 / 2 J2t=i cos ^> where 
the group velocity is w Q (k) = t/(2d) 1 / 2 sin k a . Further on we use an analytic representation ftp, the 
two-particle bubble in high dimensions where it can be represented via a double Gaussian integraE3 

/oo 
d\ 1 d\^6(X 1 lxaz- i )6{\2\mz2) 
-oo 

x exp U (\xx{zi) + X 2 x{z 2 ))-^ (A 2 + A 2 + 2\iX 2 X(q))\ (36) 

with x(z) = z — S(z) and A(q) = \ Y^l=i cos ilv), and the Heaviside step function 9(x). Only parallel 
conductivity survives on a hypercubic lattice and the integrals over the Brillouin zone of squares of the 
velocity factorize, i. e., 
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5> a (k) 2 G CT (k)G T (k + q) = ^Xar(q). 



(37) 



Using the above representations and simplifications in Eqs. ( p5| ) and (|34|) we obtain after straightfor- 
ward manipulations an explicit asymptotic formula for the conductivity in high spatial dimensions 



!_ p _ l _h^n |2\ 

e 2 Ke CTQa " 2d 



(|ImG + | 2 ) 



t 4 \ (G 2 \ 4 <" 
Re 



16rf 2 



Gl 



1 



(l + G+(£ + -^)) 5 



|l + G + (S + -l/ l )l 2 / Q , 



G\ 



- 1 



IG+I 2 



1 



|i + G + (E + -y,)l 2 / a „ 



- 1 



(1 + G+(E + -V0) 2 



(38) 



where we used abbreviations (G+) = A^ -1 G(k, Ep + i0 + ), {G\) = N^ 1 J^k G(k, -Ef + «0 + ) 2 , and 
E + = Y1(Ef +i0 + ). The one-particle propagators are the CPA ones calculated with the self-energy from 
the local approximation (|l2|). Conductivity (^8|) depends explicitly only on the disorder distribution and 
the dimensionality The vertex contribution, proportional to t^/d 2 , has negative sign and we cannot 
guarantee positivity of the conductivity. Although the vertex corrections are negligible in the limit 
d — > oo they may turn the overall sign of the conductivity negative for a fixed finite dimension d. 
Whether conductivity (|3^) for a fixed dimension is positive or negative depends on the band filling and 
the disorder distribution and strength. Vanishing of the asymptotic conductivity in this representation 
does not indicate Anderson localization but merely a limit upon the dimension below which formula d3§| ) 
cannot be applied. 



B. Mean-field expression for vertex corrections to the electrical conductivity 

The limit of high spatial dimensions is usually used in order to set a mean-field approximation for a 
physical quantity. We showed in the preceding subsection that the vertex corrections are asymptotically 
less important than the one-electron, single-bubble term. But the vertex corrections may, in finite di- 
mensions, turn static optical conductivity negative. Expression ( |3~4j ) with the asymptotic solution for the 
vertex function in high dimensions is hence unsuitable for serving as a mean-field approximation, since 
physical consistence of the result is not guaranteed. We can take the leading single-bubble term as a 
mean-field approximation for the conductivity as is actually common in the literature. Or, when we are 
interested in the impact of vertex corrections, we have to use representation ( |35| ) and evaluate its leading 
asymptotic behavior in high spatial dimensions as suggested recently.Eil 

The limit of high spatial dimensions enables one to solve the Bethe-Salpeter equation in the electron- 
hole channel explicitly. We need only its leading asymptotic order. Note that only the nonlocal part of 
the vertex A e ' l (k, k;k' — k) having odd parity with respect to reflections in k and k' contributes to the 
conductivity. We hence have to take its leading asymptotics into account. We derive it if we solve the 
Bethe-Salpeter equation in the electron-hole channel for A eh instead for T, i.e., A 6 ' 1 = T {• [GGr] + 1} . 
Using vertex T from Eq. ( p3|) one finds in the order 0(l/d) 

A eh (k 1 ,z 1 ,k 2 ,z 2 ;q) = A(z 1 ,z 2 ) 

+ (1 - A(zi, z 2 )G( Zl )G(z 2 )) 2 [r(k 1; z u k 2 , z 2 - q) - T eh (k 2 - ki; z u z 2 )} . (39) 

In the limit d — > oo the momentum convolutions decouple. To derive the leading asymptotic contribution 
from the nonlocal part of A eh to the conductivity we have to calculate the momentum convolutions at 
the level of order 0(l/<i) so that the velocities appear in squares and the momentum integrals do not 
vanish. Keeping only the leading-order terms in the expansion of the denominator in Eq. (|35| ) we end up 
with a mean-field-like expression for the dc-conductivityEil 
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where (v 2 a G a G T ) = N' 1 £ k w Q (k) 2 G CT (k)G T (k) and 



The irreducible vertex A eh is taken from Eq. (39). The one-particle propagators are those from the local, 
coherent-potential approximation, since it is the local i rredu cible vertex A that determines the leading 
high-dimensional asymptotics of the self-energy in Eq. ( |28a|) . However, only terms with odd symmetry 
with respect to time inversion contribute to the vertex corrections to the conductivity and we have to go 
to the leading order of the nonlocal part of A eh to substantiate them. If one resorts to the local vertex 
A in Eq. (^) one recovers the CPA conductivity of the single electron-hole bubble. 

Conductivity (|4^) with the CPA one-particle propagator can be called a mean-field approximation for 
the conductivity with vertex corrections, since the vertex corrections are determined from the asymptotic 
limit of high spatial dimensions. It can be applied in finite dimensions for d > 2 with the appropriate 
lattice structure determined by the Brillouin zone. Note that the vertex A eh from Eq. (p9|) contains a 



diffusion pole from the electron-electron channel and the integral in Eq. (40b) diverges in d < 2. Hence the 
mean-field description of vertex corrections break down there. In these low dimensions, where Anderson 
localization is expected, we have to use the full two-particle self-consistent parquet approximation in 
order to take properly into account the influence of the diffusion pole. 



V. BINARY ALLOY: HIGH-DIMENSIONAL APPROXIMATION 



Mean-field expression for the conductivity with vertex corrections ( p0[ ) demands the evaluation of 
momentum integrals over the Brillouin zone to determine the averaged velocity and the derivative of 
the vertex A eh . This must be done numerically for each particular lattice. To assess qualitatively the 
influence of the vertex corrections on the conductivity we can further simplify the mean-field expres- 
sion in that we approximate the momentum integrals by their high-dimensional asymptotics, i. e., we 
replace the mom entum integrals by integrals with the density of states. We use the same steps as in 
Subsection [II Cj when deriving the high-dimensional asymptotics of the conductivity. Using vertex ( J3S| ) 
we straightforwardly derive 



ReaQa -(Ld)5l + I(G CT G T >A ffT (l 



(-ar)(G a G T ) 



A aT G a G T ) [ laT {Gl){GD - 1<ja {Gl) 2 - irAGlY 



(41) 



where (G a G T ) is defined as in Eq. fl3S|). The vertices A and 7 are calculated in the local approximation 
from Subsection II A , We can explicitly sum over the indices a, r and obtain 



1^,6 CJ ' (xot 



\4ird J 



(\G+ 


to 


l + |i(|G + P)A + _(l-A + _|G + P) 


7 + -|<G2)| 2 -2Rc( 7++ (G^)" 



-Re 



(42) 



For the explicit numerical calculations we choose a binary alloy with the site-diagonal disorder distri- 
bution 



p(V) = x6(V - A) + (1 - x)S(V + A) , 



(43) 



To simplify the relation between the one and two-particle functions we further choose a model density of 
states of a d — 00 Bethe instead of a hypercubic lattice. It is characterized by an equation 



G(z) 



1 



z-G(z) 



(44) 



where we set the hopping t = 1 with the energy band E S (—2, 2). The self-energy for complex energies 
in this model is then determined from a cubic equation 
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= G(z) 3 - 2zG{zf + (1 - A 2 + z 2 )G(z) - (z - (1 - 2x)A) 



(45a) 



with 

E(z) = z - G(z) - G(z)- 1 . (45b) 

We resort to real frequencies E a = E + ia{) + and take advantage of explicit representations of two- 
particle functions via one-particle local propagators 

(G*G r ) = 1 G f T r (46) 

and 



G a G T 



(E — G a — A)(E — G T — A) - G ff + A)(£ - G T + A) 



We use the above formulas together with a numerical solution to the self-energy from Eq. (Ea) in 
Eq. (E^) to reach quantitative results for the electrical conductivity. In particular we are interested in 
the impact of vertex corrections onto the CPA conductivity. We set the formal parameter of the lattice 
dimension d — 3 in Eq. . 



x=0.5,A=0.7 x=0.5,A=1.2 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 



E/3 E/3 

FIG. 1. Real part of the static conductivity as a function of energy for concentration x = 0.5 and two values of 
disorder strength A = 0.7 and A = 1.2. 

Fig. |l| shows the CPA and the mean-field conductivity with vertex corrections as a function of energy 
for a fixed concentration x = 0.5 and two values of the variance of the random potential (disorder 
strength). We see that the vertex corrections in general lower the single-bubble conductivity. There are, 
however, situations where the vertex corrections may lead to an increase in the conductivity. It is mostly 
due to proximity of nonanalyticities in the local CPA vertex. The real part of A ++ for A = l/v2 goes 
through zero and displays a pole for x = 0.5. This singular behavior causes sharp mobility peaks at the 
band edges and irregularities near half-filling, E — 0. Such a behavior is observed only for A w 1/V2 
and x m 0.5. Fig. |^ shows the same for an asymmetric concentration x = 0.3. The mobility edges are 
now less pronounced and the conductivity fluctuations inside the energy band are no longer symmetric. 
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x=0.3, A=1.2 
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FIG. 2. Real part of the static conductivity as a function of energy for concentration x = 0.3 and two values of 
disorder strength A = 0.7 and A = 1.2. 

Due to the structureless density of states of the Bethe lattice in d = oo, vertex corrections do not alter 
the CPA conductivity significantly apart from the special situations influenced by the singularity in the 
CPA vertex A ++ . Fig. || shows the conductivity for the half- filled band as a function of concentration x. 
There is almost no significant difference between the conductivity with and without vertex corrections, 
except for concentrations x « 0.5. An analogous picture is obtained for the conductivity as a function of 
the disorder strength A, where the differences are significant only around A w l/\/2, Fig. ||. 
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FIG. 3. Real part of the static conductivity as a function of concentration for the half-filled band and A = 0.7. 
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FIG. 4. Real part of the static conductivity as a function of disorder strength for the half-filled band and 
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VI. CONCLUSIONS 



We have developed a diagrammatic theory for constructing systematic approximations to nonlocal 
two-particle vertex functions of noninteracting electrons moving in a random potential. The underlying 
idea of our approach is to treat separately diagonal, local and off-diagonal, nonlocal elements of the 
two-particle vertex. To this purpose we used the asymptotic limit of high spatial dimensions. In the 
strict d — oo limit only the local one-particle propagator is relevant and the solution contains only 
single-site scatterings and reduces to the coherent-potential approximation. Beyond this limit we utilized 
ambiguity in the definition of the two-particle irreducibility. We classified the nonlocal contributions 
to the two-particle vertex according to the type of the two-particle irreducibility they belong to. Alike 
many-body theories there are three topologically inequivalent irreducibility channels according to what 
pairs of propagators interconnect spatially distinct two-particle scatterings on the random potential. 
Representing the two-particle vertex via Bethe-Salpeter equations and irreducible vertices in each channel 
and utilizing the topological inequivalence of these representations we derived a closed set of coupled 
(parquet) equations for the irreducible vertices. The irreducible vertices from the parquet equations were 
used in an integral form of the Ward identity to determine the self-energy of the parquet solution. In this 
way we completed the parquet equations to an approximation consistently determining all one and two- 
particle functions. The input to the parquet approximation are the local CPA one-particle propagator 
and two-particle vertex. Neither the form of the parquet equations nor the Ward identity may cause 
unphysical nonanalyticities in the solution. Solutions to the parquet approximation hence inherit the 
analytic properties of the CPA input and are free of spurious, unphysical behavior. 

The proposed diagrammatic implementation of nonlocal corrections to the CPA aims primarily at 
improving the CPA two-particle functions on a long-range scale. Although there is no obvious small 
parameter controlling the nonlocal corrections to the CPA, systematic improvements of the local approx- 
imation are controlled via diagrams to the completely two-particle irreducible vertex /. In the parquet 
approximation the input I = -f, the local vertex from Eq. (|lq). A first correction AI to the input of the 
parquet approximation is proportional to 7 4 (G°^) 6 , where is the off-diagonal element of the CPA 

propagator, see Fig. ||. It means that in the weak-disorder limit the parquet approximation is exact up to 
V 7 whereas CPA only to V 3 in powers of the random potential. The parquet approximation represents a 
significant systematic improvement of the CPA in the weak-disorder limit. We showed that beyond the 
weak-disorder limit the parquet approximation contains the exact asymptotics of the two-particle func- 
tions up to dr 2 . This fact we used in proposing a mean- field approximation for the electrical conductivity 
with vertex correction. 



ak 



ak + q 




a'k' + q 
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>i 
>> 
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FIG. 5. Lowest-order correction to the completely two-particle irreducible vertex I aa i = 7 aa / from the parquet 
approximation. The double-dashed line is the CPA local vertex r y aa > . The internal fermion lines stand for the 
CPA off-diagonal propagator . 



It is not, however, the weak-scattering limit where the assets of the parquet approach to disordered 
systems lie. The parquet approach is in particular appropriate for apprehending spatial quantum coher- 
ence and backscattering effects. The parquet approximation even in the first iteration (high-dimensional 
asymptotics) contains infinite number of "crossed" two-particle diagrams. The weak localization with the 
diffusion pole and a long-range spatial coherence from backscatterings are included. Most importantly 
however, the parquet equations for the irreducible two-particle vertex functions are fully self-consistent. 
They can adequately deal with poles and divergences in the vertex functions and can hence significantly 
change the diffusive character of the electronic transport in low dimensions and strong disorder. 

We hope that the presented parquet approach can bridge the gap between the mean-field coherent- 
potential approximation on the one side and localization theories on the other side. Whether the parquet 
equations can actually describe the localization transition must be decided by solving the full set of self- 
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consistent coupled integral equations for the vertex functions and the self-energy. This has been left for 
future research. 
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APPENDIX A: ALGEBRAIC DERIVATION OF WARD IDENTITIES 



We use algebraic identities to prove equation (|26|) and the special case of fl27|) with ki = k 2 . 

The two-particle function can be defined as a matrix element of a tensor (direct) product 



( 2 ) ^ ^ ^ 

G\- kl (zi,z 2 ) — {{ik\G(zi)®G(z2)\jl))av where the resolvent operator is defined as G(z) = 



zl-t-V 



and the basis vectors are Wannier orbitals at the lattice sites. Using this lattice-space representation we 
define a "projection" onto the one-particle subspace by equaling the basis states from the left and right 
Hilbert space. When we sum over one set of indices (k = j) we reduce the direct product from the 
two-particle Hilbert space to an operator multiplication in the one-particle Hilbert space. Ward identity 
(|2q) is then a consequence of an operator identity 



zxl - H 



Zol - H 



= \Z2 - Z\\ 



1 r ^ ~n -1 

zil-H - z 2 l-H 



(Al) 



where H = t + V . Identity (Al) holds for each configuration of the random potential V and after its 
configurational averaging we obtain Eq. (p6|). Ward identity (26) is hence a consequence of completeness 
of the eigenstates of the Hamilton operator. Completeness of the eigenstatcs reflects conservation of 
probability and is a necessary prerequisite for conservation of energy and other physical quantities. 

To prove the Ward identity (^) with ki = k2 we use the Bethe-Salpeter equation in the electron-hole 
channel 



Gfl kl {z l ,Z2) = G ij (z 1 )G u (z2)+ £ G H/ (z 1 )G w (z2)At jl , mi (z 1 ,z 2 )Gf> M ,(z u z 2 ). 

i'j'k'l' 

We multiply it with the inverse one-particle propagators from left and right and obtain 

£G i :, 1 (z 1 )G H , 1 (z 2 )Gg fci ,(z 1 ,2 2 ) = SyS* + £ ^^(zi, z 2 )Gf jMt ( Zl , z 2 ). 

i'V j'k' 

Summing over the intermediate indices and using Eq. ( |2^ ) we find 

1 



J2 G u'^)GuH^)Gg jV (zuz 2 ) 

i'jl' 



Az 



[G a \zi)-G-\z 2 )] 



We insert Eqs. (A4) and (|26|) in Eq. (|A3[) and obtain the desired identity for the self-energy 



It reads in momentum space 

E(k, zi) - S(k, z 2 ) = 1 Aeh (^ ^; k' - k) [G(k', zi) - G(k', z 2 )] 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



Analogously we can derive a Ward identity for the vertex function from the electron-electron channel 
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S(k, zi) - E(-k, z 2 ) = 1 J] A ee (k, * x> k', z 2 ; k' - k) [G(k', z x ) - G(-k', z 2 )] . (A7) 



TV 
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